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Abstract 

We reproduce the mid-infrared to radio galaxy counts with a new empirical model based on our current 
understanding of the evolution of main-sequence (MS) and starburst (SB) galaxies. We rely on a simple Spectral 
Energy Distribution (SED) library based on Herschel observations: a single SED for the MS and another one 
for SB, getting warmer with redshift. Our model is able to reproduce recent measurements of galaxy counts 
performed with Herschel, including counts per redshift slice. This agreement demonstrates the power of our 2 
Star-Formation Modes (2SFM) decomposition for describing the statistical properties of infrared sources and 
their evolution with cosmic time. We discuss the relative contribution of MS and SB galaxies to the number 
counts at various wavelengths and flux densities. We also show that MS galaxies are responsible for a bump in 
the 1.4 GHz radio counts around 50 juJy. Material of the model (predictions, SED library, mock catalogs...) is 
available onlin^J. 

Subject headings: galaxies: statistics — galaxies: evolution — galaxies: star formation — infrared: galaxies 
— submillimeter: galaxies 



1. INTRODUCTION 

Recent observational studies have shown that two distinct 
star-forming (SF) mechanisms are required to describe the SF 
galaxy population. The so-called SF main sequence (MS) is 
composed of secularly-evolving galaxies that display a tight 
correlation between stellar ma ss (M+) and star formation rate 
(SFR) at a given redshif t (e.g. lElbaz et al.ll2007HNoeske et al l 
120071: iDaddi et al.ll2007l) . This population accounts for -85% 
of the star formation rate density (SFRD) i n the Universe 
dRodighiero et al.1 1201 U ISargent et all 120121) at z<2. The 
rest of the star-formation budget is provided by starbursts 
(SB), i.e. galaxies with very high specific star formation 
rates (sSFR=SFR/M A ) probably induced by re cent mergers 
(e.g. lElbazetaLlboi U iRodighiero et all 12011b . Recently, 
Sargent et al. (2012, S12 hereafter) showed that infrared 
(IR) luminosity functions (LF) can be reproduced by jointly 
considering the mass function of SF galaxies (SFMF), the 
evolution of the sSFR of MS galaxies, and its distribution 
at fixed M*, with a separate contribution from MS and SB 
galaxies. 



1 Laboratoire AIM-Paris-Saclay, CEA/DSM/M'u - CNRS - Universite 
Paris Diderot, CEA-Saclay, Orme des Merisiers, F-91191 Gif-sur-Yvette, 
France, Email: matthieu. bethermin@cea . fr 

2 Department of Physics, University of Oxford, Keble Road, Oxford 
OX1 3RH, UK 

3 Department of Physics, McGill University, 3600 Rue University, Mon- 
treal, Quebec H3A 2T8, Canada 

4 UPMC Univ. Paris 06, UMR7095, Institut d'Astrophysique de Paris, 
75014 Paris, France 

5 CNRS, UMR7095, Institut d'Astrophysique de Paris, 75014 Paris, 
France 

6 Laboratoire d'Astrophysique de Marseille, OAMP, Universite Aix- 
Marseille, CNRS, 38 rue Frederic Joliot-Curie, 13388 Marseille Cedex 13, 
France 

7 Department of Physics & Institute of Theoretical and Computation 
Physics, University of Crete, 71003 Heraklion, Greece 

8 Institut d'Astrophysique Spatiale (IAS), batiment 121, Universite 
Paris-Sud 11 and CNRS (UMR 8617), 91405 Orsay, France 

9 Department of Physics & Astronomy, University of British Columbia, 
6224 Agricultural Road, Vancouver, BC V6T 1Z1, Canada 

' at http://irfu-i.cea.fr/Pisp/matthieu.bethermin/ 



Wavelength-dependent galaxy number counts are an 
additional, important constraint for evolutionary models 
of in frared galaxies. While purely semi-analy tical models 
(e.g. lLacev et all 120 lOt ISomerville et al.l 120 12f) struggle to 
reproduce infrared (I R) number counts, phenomenological or 



hybrid models (e.g. Bethermin et al. 
1201 U iRahmati & van der Werfl 1201 1 



lLapi et all 1201 lh fare 



et a. 
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better but are in general descriptive and use an evolution of 
the luminosity function which is not motivated by physical 
principles. However, these recent models which reproduce 
the total counts passably, are excluded at >3 cr by the 
recent Herschel measurements of counts per redshift slice 
dBerta etalJl20lU iBethermin et alJl2012bl) . This shows how 
important redshift-dependent constraints are to accurately 
model the evolution of galaxies, and motivates the develop- 
ment of a new generation of models. 

We present a new model of IR galaxy counts which builds 
on the 2-Star- Formation-Mode framework (2SFM) S12 
introduced. This fiducial model is intuitive and based on 
our current observational knowledge of the evolution of 
MS and SB galaxies. All model parameters are constrained 
by external datasets and require no additional fine-tuning. 
We assume a Salpeter initial mass function and a WMAP-1 
cosmology. 



2. MAIN INGREDIENTS 

Our model is based on four main ingredients, which are 
sufficient to reach a good agreement with IR source counts 
(see Sect.|4]and grey line Fig. [3]): 

• evolution of the MS with redshift, 

• decomposition of the sSFR distribution at fixed M* into 
MS and SB modes, 

• evolution of the SFMF with redshift, 
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• Spectral Energy Distribution (SED) libraries for MS 
and SB galaxies. 

Additional ingredients, which are of lesser importance, are 
presented in Sect. [3] 

2.1. SFR distribution 

A key ingredient of the S12 approach is the probabil- 
ity distribution of sSFR at fi xed M» for SF galaxies based 
on observations presented in Rodighie ro et alj d201 ll) . It is 
parametrized as a double log-normal decomposition of MS 
and SB: 



p(log(sSFR)) oc exp 



( (log(sSFR)-log(sSFR MS )) 



2n 



2 <s 



+r SB x exp 



( Qog(sSFR) - log(sSFR M s) - B SB )) 
2cr L 



2\ 



(1) 



where cr MS and cr s B are the dispersion in the sSFR of the MS 
and the SB populations. B$ b is the average sSFR-boost for SB 
galaxies. We assume that these three parameters do not evolve 
with M+ and redshift, as suggested by S12 who reproduce the 
z~0 IR LF under these assumptions and with the distribution 
calibrated at z~2 (see Table Q]for parameter values adopted). 
sSFR MS varies with M* and redshift according to 



sSFR MS (z,M*) =sSFR MlS ,oX 



ft MS 



(2) 



x(l +min(z,z evo )) 



7ms 



where sSFR M s,o is the sSFR at z=0 for M* = 10 11 M Q and 
B MS parametrizes the dependance of sSFR on M*. y MS de- 
scribes the evolution of the normalization of the MS out to 
redshift z evo = 2.5 where this evolution fl attens according 
to observations (e.g. Gonza lez et al.ll20loh . The values of 
these parameters, chosen based on measurements summarized 
in Fig. [The, are listed in TableQ] S12 also present evidence 
for a weak redshift evolution of tsb, the relative amplitude of 
SB sSFR log-normal distribution compared to MS one (or, 
equivalently, of the relative SB-contribution to the S FRD, 
see Fi g. [pi), in agreement with the model of Hopki ns et al.l 
(2010). Here we define the redshift evolution of r SB as: 

fssiz) = r SB ,o x (1 + mm(z,ZsB)) ysB , where z S b = 1, (3) 

in order to broadly reproduce the trends suggested by these 
two studies (see Fig. [TJl). The impact of this evolving r$s 
is negligible, barring a ~20% decrease of 70yum counts 
compared to a constant r$B- 

Another important ingredient of our model is the evolution 
of the SFMF Observations are well-described by a Schechter 
function 

dN (MA- am I MA 

d> = =d>h(z)x\ — xexp x ln(10) 

V dXogiM*) m; \M b J P \ Mb) M h 

(4) 

with a redshift-invariant cha racteristic mass M /, and faint-end 
slope auF, in keeping with iPenjTet al. (2010). <pb, the char- 
acteristic density, is constant between z=0 and z=l but de- 
creases at z>l as 
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Fig. 1 . — Redshift evolution of selected model parameters and derived quan- 
tities; model conventions are represented using a solid line; (a) - density at 
the break of the mass function; (b) - sSFR at M* = 5 X 10 9 M Q . (c) - sSFR at 
M t = 5 X 10 10 M Q (data in a, b, and c from a compilation of Sargent et al. in 
prep.); (d) - contribution o f SB to the total star formation rate density (data 
from S12; grey region from lHopkins et alJ20iq) . (e) - SFRD . The MS (SB) 
contribution is represented by a dotted (dashed) line (data: Bouw ens et alj 
<2007h .lRodighiero et al. 2011, Magnelli et al. 2011 and Kari m etaU2011l) . 



\og((p b ) = log(0&)(z < 1) + Tsfmf(1 - z). 



(5) 
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TABLE 1 

Summary of parameters of our fiducial model. 



Parameter 


Description 


Value 


Reference 


Distribution of sSFR 


B SB 


Boost of specific star formation rate in SB (in dex) 


0.6 


Sargent et al. (2012, hereafter S12) 


cms 


Width of the MS log-normal distribution (in dex) 


0.15 


S12 value minus 0.05 dex for artificial scatter (Salmi et al. 2012) 




Width of the SB log-normal distribution (in dex) 


0.20 


S12 value minus 0.05 dex for artificial scatter (Salmi et al. 2012) 




Evolution of the 


main MS 




sSFR M5 - 


sSFR on the MS at z = and M* = 10 11 Af (in log(yr')) 


-10.2 


S12 based on a fit of literature data 




Slope of the sSFR-M* relation at a given redshift 


-0.2 


Rodiehiero et al. (2011) and S 12 


7ms 


Evolution of the normalization of the MS with redshift 


3 


compilation of measurements of sSFR (see Fig.fT} 


Zevo 


Redshift where the MS normalization stops to evolve 


2.5 


compilation of measurements of sSFR (see Fig.fT} 


Evolution of the fraction of SB 


r SB,0 


Relative amplitude of SB log-normal distribution compared to MS 


0.012 


S 12 and Hopkins et al. (2010) (see Fig.fTl 


ysB 


Evolution of starburst fraction with redshift 


1 


S 12 and Hopkins et al. (2010) (see Fig.fT} 


ZSB 


Redshift where the starburst fraction stops to evolve 


1 


S12 and Hopkins et al. (2010) (see Fig.fT} 


Star-forming mass function (SFMF) and its evolution 


M b 


Stellar mass at the break of the SFMF (in log(M Q )) 


11.20 


Ilbert et al. (2010) and Peng et al. (2010) 


a 


Faint-end slope of the SFMF 


1.3 


Ilbert et al. (2010) and Pens et al. (2010) 


<fa(z < l) 


Number density at the break of the SFMF at z < 1 (in log(Mpc~ 3 )) 


-3.02 


Sargent et al. (2012) (see Fig.m 


TSFMF 


Evolution of the density of SFMF at z > 1 


0.4 


Extended from Sargent et al. (2012) (see Fig.fT} 


Evolution of SEDs 


(U)mso 


Mean radiation field in local MS galaxies 


4 


Magdis et al. (2012) and Fig.ff] 


7U.MS 


Evolution of the radiation field in MS with redshift 


1.3 


Magdis et al. (2012) and Fis.l2l 


Z(U),MS 


Redshift where ((/) in MS flattens 


2 


Magdis et al. (2012) and Fig.ff] 


(U)SB,0 


Mean radiation field in local SB galaxies 


35 


Magdis et al. (2012) and Fig.ff] 


7U.SB 


Evolution of the radiation field in SB with redshift 


0.4 


Magdis et al. (2012) and Fig.ff] 


Z(U),SB 


Redshift where (U) in SB flattens 


3.1 


Magdis et al. (2012) and Fig.ff] 


Contribution of AGNs 


A-AGN 


Normalization of the probability distribution of AGN contribution 


0.0025 


Airdetal. (2012) 


Pagn 


Slope of the probability distribution of AGN contribution 


-0.7 


See Sect.fJ 



The fiducial values (chosen from Fig. [T^) of the MF-related 
parameters are also listed in Table Q] 

The star formation history implied by our evolutionary 
formalism is shown in Fig. QJ. The star formation rate 
density increases from z=0 to 1, flattens between z — 1 
and z = z evo = 2.5 and decreases with redsh ift at z > z evn , 
match ing the infrared measurem ents of Magnel li et alJ 
(1201 ll) and Rodighie ro et alJ (1201 lh . the radio measurements 
of iKarim et al.l d201 ll) . and the optical measurements of 
Bouwens et al. (2007) at high redshift. The SFMF is quite 
uncertain at z > 4, but this has little impact on the counts. 
In our model, the SFRD is dominated by MS galaxies at all 
redshifts. 



2.2. SEDs 

We use a ch aracteristic IR SEP template for MS and SB 
based on fits of Drai ne & Lil ((2007) mode ls to Herschel obser- 
vations of distant galaxies as presented in Magdis et al.l d20 1 21 
hereaf ter M12). While conceptually similar to Elbaz et al. 
( 1201 ll) . however, our templates evolve with redshift following 
the finding of M12 that the mean radiation field (U) (which 
correlates with dust temperature) is more intense at high red- 
shift: 

(U) = {U) x(l+rmn(z,z {U) )y u . (6) 

Here (£/)o is the mean radiation field in local MS galaxies, ju 
a parameter determining its evolution with redshift, and Z(u) 
the redshift where (£/) flattens. This evolution is different in 
MS and SB galaxies (see Fig. [2] and TablefTJ. This evolution 
is caused by the evolution of SF efficiency and metallicity 
with redshift (Ml 2), and is required to reproduce source 



counts. For example, if we used the z=l (z=0) MS template 
for all redshifts, we would overestimate the counts by about a 
factor of 2 (1) at 70 /im and 2 (10) at 1 . 1 mm. For re f erence , 
if we use the MS and SB templates of lElbaz et al] (12011b . 
we overpredict the mm counts by a factor of 10 at all fluxes 
and underpredict the 100 /vm counts by ~30%. To reproduce 
the 24yum counts, it is crucial to use distinct SB templates 
with less mid-IR emission than in MS galaxies. The SEDs of 
MS and SB galaxies used in our model are shown in Fig. [2] 
We introduce a relative dispersion on (U) of 0.2 dex for both 
MS and SB (Ml 2), which has little impact on the counts 
(<10%), except in the millimeter domain (+20%). In this 
approach, the increasing mean dust temperature with infrared 
luminosity (L/r) at a given redshift is caused by a higher 
fraction of SB galaxies at higher L/«. 



3. REFINEMENTS 

3.1. Dust attenuation 

To reproduce IR number counts we have to link SFR and 
L/r. For obscured SF galaxies, the bulk of the UV light emit- 
ted by young stars is absorbed by dust and re-emitted in th e 
IR (SFR/r/L/r = K = 1.7xlO- 10 Mnvr 1 Lp 1 . rKennicuttll998l) . 
In less massive galaxies, the attenuation is smaller and a sig- 
nificant part of the SF can be detected in UV. The total star 
formation can then be divided into an uncorrected UV and an 
IR component (S FR - S FRuv + S FRjr). The mean ratio 
between these two comp onents, ri sno, varies with M+. Here, 
we apply the relation of Pan nella et al.l (120091) : 

( SFR, R \ (MA 

^- 2 - 5i °4sF^) =4 -° 7xiog y- 39 - 32 ' (?) 



4 



Bethermin et al. 



1.00 



0.10 




0.01 



1.00 



0.10 



0.01 



Main sequence 



z=0.0 (<U>=4.0) 
z=0.5 (<U>=6.8) 

z=1.0 (<U>=9.8) 
z=1.5 (<U> = 13.2) 
z=2.0 (<U> = 16.7) 




Starburst 




10 100 

Wavelength [urn] 



1000 



100 



A 
V 



10 



I I I 

- OMS 
: A SB 








i i i i 



Fig. 2. — SEDs (normalized to have Ljr = 1 L Q ) for MS (top) and SB (mid- 
dle) used in our model. Evolution of (U) parameter with redshift in MS (blue) 
and SB (red) galaxies (data from M12). 



and as sume redshift-invariance, as suggested by Sobr al et all 
d2012l) and Pannella et al. (in prep.). The IR luminosity of the 
galaxies, L/r is thus given by: 



SFR 



■IR 



SFR 



JQ0.4XH500 



SFR 



K 



K 1 + 10° 



4x)-i . 



K 



X/i 



where /^ F (M*) 



(8) 



S FRjr/S FR goes to at low mass and 
1 at high mass. This correction implies a flatter IR LF at 
the faint end as compared to the SFMF at the low-mass end 
and prevents an excess in the counts at faint flux densities. 
Although a small part of the IR emission is due to dust heated 
by old stars, especially at low-z, we consistently reproduce 
z = 0-2IRLF(S12). 



3.2. AGN contribution 

Active-Galactic-Nucleus (AGN) activity is potentially im- 
portant when modeling mid-IR counts. We statistically asso- 
ciate an AGN contrib ution, represented by th e average intrin- 
sic SED templa te of iMullanev et al.l (1201 ll) . to each galaxy 
based on its Ljr. lAird et a 1J (I2012D showed that the Eddington 
ratio rEdd (bolometric luminosity L? , N over Eddington lumi- 
nosity) of AGN at z<l follows a power-law probability dis- 
tribution function (PDF) with redshift-dependent nor maliza- 
tion. Based on the results of Mullane y et alJ (120121) - who 
report a coincident cosmological evolution of the averages of 
specific black hole (BH) growth (Mbh/Mbh, where Mbh is 
BH mass) and sSFR over 0.5<z<2.5, a fact that impli e s con- 
stant Mbh/M+ ratios - , we can express the lAird et al.l d2012l) 
results in terms of a distribution of ratios of bolometric lumi- 
nosities (r AGN = Lfg N IL S j[) from AGN and SF with redshift- 
independent normalization: 



p{rEdd) = C(z) x r\ 
where we recall that 



AGN 

Edd 



p(r AC N) = A ACN x rf™, (9) 



M, 



BII 



tAGN 

Hoi 



jAGN 



M, 



BII 



j AGN 



rSF 
^IR 



(10) 



The last ster 



uses the M*-SFR correlation. f3 ACN = -0.7 
agn is based on the normal- 



tst step uses tne M^-a rK 

comes from Uirdet al.1 (120121) A ,„ 

ization of the lAird et al.l d20 12b relation and includes a scaling 
factor for the conversion between r^dd and r A cN PDFs. This 
scaling relation assumes a mean ratio between black hole and 
stellar mass of 0.0015 (Mullanev et all I20l2|). plus a mean 
ratio b e tween Lfj™ and L j™ c alibrated fromlMullanev et alJ 
(l20Tlh . iLutz et alJ d2004)°and IVasudevan & Fabianl (120071 
In order to normalize this PDF, we place a cut at Aem = 1 
and choose a lower bound such that f p(r A cN)dr A cN = 1 ■ We 
emphasize that Eq.|9]implies a correlation between AGN and 
SF activity only in an average sense, while preserving a large 
dispersion for individual objects consistent with observations. 
Full details of our AGN-treatment will be presented in a 
future paper. The AGN contribution is significant (> 10%) 
only at 24fj.m above 3mJy (see Fig. [3]) and negligible at 
longer wavelengths (< 2%). 



3.3. Magnification caused by strong lensing 

Having computed the IR LF, split into MS and SB contri- 
bution as in S12, we include the effect of the strong (fi > 2) 
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lensing dNegrello et alj|2007ll2010t) on these two LFs: 
d 2 N , r dPQi,z) d 2 N 



dlogL [R dV' 



f 



dlogjj dlogL IR dV 



! initial 



dlogfx, 



(11) 



where p. is the magnification, the magnification PDF in 
the iHezaveh & Holder! (1201 lh model, and 



-t, — ; — w-. the lumi- 
nosity function. These lensed sources contribute ~2Q% to 
(sub-)mm counts around 100 mJy. 

4. RESULTS 
Number counts are computed according to: 

d 2 N r r 

—-(S,z,A)= V dr ACN d(U) 

dSdZ type=1^S.SB| J <^ J ^ 

d 2 N, me 

x^^(z,L IR (S, type, {U),z}) 



dLj R dV 
dV 



s£ W fc X) + r AC N x S£™{z, A) dz 



p(fAGN) p«t/)k,type) 
(12) 



Here S 



type,(£/> 



is the flux of a L IR — 1 L source of a given 



type (MS or SB) and a given (U) in a given filter. S^^Z is 
the same quantity, computed using the Mull anev et al.l (201 1) 
AGN template. Note that the filter shape is taken into account 
for the calculation of S v . /?(?"agn) is provided by Eq. [9] and 
p({U)) is a log-normal distribution with a width of 0.2 dex 
(see Sect.EJ. 

We compare the predictions of our model with measure- 
ments of differential galaxy counts from 24 pm to 1.1mm 
(see Fig.0. Spitzer and Herschel counts are well reproduced, 
showing the effectiveness of our new approach. Note, 
however, a 10-20% (~ 2cr) excess at 24 /mi between 400 fiiy 
and 2mJy, and a ~20% (~ 2cr) excess at the faint-end at 
70 and 160//m (< 1 mJy and < 5mJy, respectively). The 
BLAST and SPIRE counts at 250, 350, and 500 pm are 
globally well reproduced. Neverth eless, the model slightly 
overpredicts the three last points of Bether min et al.1 ([2012b) 
(in red). As discussed by these authors, this could be related 
to an under-density in GOODS-N. The contribution of 
lensed sources broadly agre es with the measurements of 
iGonzalez-Nuevo et ail (|2012) at 350/rni. At 1.1mm, our 
model nicely agr ees with the combined number counts of 
IScott et all (l2012 i frcL except for the faintest point, originating 
from 1-2 cr sources and potentially poorly de-biased. 

S ince the LF evolution an d number counts may be degener- 
ate dBethermin et al.ll2012bl) . galaxy counts split per redshift 
provide a powerful test of the validity of our model (note that 
S12 demonstrated that bolometric IR LF are reproduced at 
z<2.5). This observable is close to the monochromatic LF, 
but requires fewer corrections (K-corrections, V max ) which 
could bias the results (possible biases from photometric 
redsh ifts and source identification are discussed in Ber ta et aT] 
1201 ll and [Bethermin et al. 2012b). The comparison between 
our model and observations (Fig.[4j reveals a good overall 
agreement between predictions and data. However, we 

10 The counts showed in Fig.[3]are corrected for the bias found in their 
simulation. 



slightly over-predict the counts around 500 p]y between 
z=0.5 and z=2 at 24 pm. It could be due to a slight excess 
of PAH features around 15 pm in the SB templates. We also 
underpredict the counts at 1 00 and 1 60 pm by 1 -2 cr, probably 
due to a slight lack of warm dust in the SED templates. 
Finally, our model overpredicts by ~ 3 cr the z > 2 counts in 
the 2-6 mJy range. As explained in the previous paragraph, 
this could be due to cosmic variance, as these points rely 
exclusively on GOODS-N. 

By distinguishing between MS and SB activity, the 2SFM 
framework allows us to explore selection biases toward MS or 
SB objects in surveys probing various wavelengths and flux 
density regimes. MS galaxies (dotted line in Fig.0 dominate 
the number counts at all flux densities and all wavelengths. 
However, the relative contribution of SBs varies a lot with 
flux density and wavelength and is important (~30%) around 
30 mJy at 70 pm and 50 mJy at 350 and 500 pm. The relative 
contribution of SB is very sensitive to the evolution of their 
SED, which is few constrained. If (U) did not evolve with 
redshift, SBs would dominate around lOOmJy at 350 and 
500 pm and at flux densities larger than 8 mJy at 1 . 1 mm. 

Finally, by assuming a non-evolving IR-radio correlation 



(gr/* = log( , 7S £L w x%^) = 2.64) out to high redshift 
(e.g. ISargent et al.l 12010) and a synchrotron spectral slope 
a=0.8 (S v oc v~ a ), we also investigate the contribution of SF 
galaxies to radio source counts at 1.4 GHz (see Fig. [3). We 
combined our model for star-forming o bjects with the model 
of AGN-driven radio sources of Massar di et al.l d2010b . The 
result agrees with the compilation of Vemstrom et al.l (1201 ll) 
(see Fig.[3]l. According to our model, the 1.4 GHz counts are 
dominated by SF obje cts below 200 pjy, in agreement with 
the observations of e.g. [Sevmour et al. ( 2008). We predict the 
presence of a bump in the Euclidian-normalized radio counts 
around 40 //Jy which is essentially due to MS galaxies. 



5. CONCLUSION 

Our model based on the main assumption of two SF modes 
(MS and SB) is able to accurately reproduce the emission of 
galaxies integrated over most of the Hubble time as probed 
by galaxy counts from the mid-IR to radio wavelengths. 
This model contains two main ingredients: the evolution of 
MS and SB galaxies based on the S12 formalism and a new 
library of MS and SB SEDs derived from Herschel observa- 
tions (M12). Despite its simplicity, our model provides one 
of the best fits achieved so far to the number counts, including 
counts per redshift slice in the SPIRE bands, which were 
poorly reproduced by the previous generation of models. All 
these results were obtained without any arbitrary tuning of 
parameters that are not constrained by observations, contrary 
to most previous models. The decomposition into 2 modes 
of SF (2SFM), i.e. MS and SB, associated with two different 
families of SEDs, is thus a very powerful framework to 
statistically describe the dust emission of galaxies across 
cosmic time. In addition, we present a new stochastic AGN 
treatment, and also found that MS galaxies are responsible 
for a bump in the 1 .4 GHz radio counts around 50 pJy. 

This model can be combined with halo mo dels assuming a 
link be t ween SFR, M*, an d halo mass (e.g. iBethermin et alj 
l2012at I Wang et all 120121) to interpret the clustering of 
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Fig. 3. — Number counts from 24 pm to 1.4 GHz. Solid line- total counts predicted by the model; grey line - counts predicted by the simplified model (without 
refinements discussed in Sect. [5J; dotted line - MS contribution; short-dashed line- SB contribution; dot-dashed line - lensed so urces; triple-dot-dashe d line- 
difference between counts with and without AGN contribution. At 1 .4 GH z, we also pl ot the model of AGN-driven radio sources of Massardi et al. 1 2010) (long- 
dashed line) and combine it with our model of SF galaxies. Data points - Bethermin et al. (2010a) (red points at 24, 70, and 160pm), Berta et al. (201 1) (blue 
points 70, 100, and 160/im), Bethermin et al. (2012b) ( red po i nts at 250, 350 and 500um). IOliver et alj<201fl) (blue points at 250, 350 and 500um). IGienn et all 
(20101) (gTeen points at 250, 350 and 500 om), IClements et al l J2010h (yellow points at 250, 350 and 500pm),[Bethermin et al. 1 2010b) (purple points at 250, 350 
and 500pm), Scott et al. 12012) (red points at 1.1 mm) , and Vernstrom et al. (201 1) (compilation of 1.4GHz radio counts). Black dots - contribution of lensed 
galaxies at 350/im measured by Gonzalez-Nuevo et al. 12012). 

infrared galaxies and the fluctuation of the cosm ic infrared 
background (e.g. iPlanck c ollabora tion et al.l 1201 ll) . Finally, 
this model and its future extensions will provide predictions 
for the next generation of IR, mm, and radio surveys, and, 

in particular to anticipate which galaxy populations will be We acknowledge Kimberley Scott, Herve Aussel, Emeric 

preferentially detected, depending on the survey strategy Le Floc'h, Benjamin Magnelli, the anonymous referee, ERC- 
adopted. ' ' StG UPGAL 240039, and ANR-08-JCJC-0008. 
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